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Abstract 

A  method  for  obtaining  the  complex  roots  of 
polynomials  of  degree  less  than  or  equal  to  twenty 
with  real  coefficients  is  described.   Both  the  mathe« 
matlcal  and  Unlvac  aspects  are  discussed. 
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Section  lo   Introduction 

A  procedure   for  obtaining  the  real  and  complex 
roots  of  polynomials  of  degree  twenty  or  less  with  real 
coefficients  has  been  tested  on  the  Unlvac,   In  general 
the  program  produces  answers  correct  to  six  or  seven 
significant  figures  (on  the  basis  of  floating  point 
operations  with  nine-digit  numbers).   Some  outside  inter- 
vention is  necessary  when  there  are  roots  which  are  "very 
close"  to  each  other,  or  when  the  coefficients  of  the 
polynomial  have  "too  wide"  a  spread  (see  Section  3  of  this 
report) • 

The  procedure  makes  repeated  use  of  the  Euclidean 
algorithm  for  finding  the  greatest  common  divisor  of  two 
polynomials.   The  fact  that  the  Cauchy  formula  which  gives 
the  number  of  zeros  to  the  "left"  of  a  straight  line  can 
be  evaluated  by  performing  the  Euclidean  algorithm  is  de- 
rived in  Section  2, 

It  was  found  to  be  a  matter  of  mathematical  ne- 
cessity to  first  eliminate  the  occurrence  of  multiple 


1,   The  method  used  has  been  described  in  the 
literature  by  T,  C,  Fry  (Quarterly  of  Applied  Mathematics, 
Vol.  Ill,  No.  2,  July  19i|.5)  who  states  that  G.  R ,.  Stibitz 
has  programmed  such  a  procedure  for  polynomials  of  degree 
ten  or  less.   The  method  may  ultimately  be  traced  back  to 
Routh,, 
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roots.   This  is  done  by  factoring  out  the  greatest  common 
divisor  of  the  polynomial  and  its  derivative.   (Since  the 
principal  tool  used  in  the  location  of  the  roots  is  based 
on  the  Euclidean  algorithm  for  finding  the  greatest 
common  divisor  of  two  polynomials,  this  is  a  simple  pre- 
caution, ) 

The  next  step  is  to  determine  a  circle  about  the 
origin  in  which  the  roots  lie.   A  simple  estimate  of  the 
form 


max 
zl  <  


^i 


+  1  +  10'^  =  B 


^nl 


is  used,  where 

P„(z)  =  a^z^  +a„,z"  +»»»  +  a,  z  +  a^. 
n      n      n-l  1     0 

The  abscissae  of  the  roots  are  now  located  by 
successive  subdivision  as  follows.   A  record  of  the  lower 
and  upper  bounds  of  the  abscissae  of  each  of  the  n  roots 
Zt»  Zo*  •••»  ^n   ^^    kept;  for  example,  at  the  start 

-  B  <  RqZ^  <  +  B  for  i  =  1,  2,  ...,  n. 
If  the  upper  bound  G-,  and  the  lower  bound  g^  of  the 
abscissa  of  z,  satisfy  the  inequality 

%        -.^=^-   <  10-8 
|G+g|+l 


'k' 


then  the    avorcige  ••7   (g^-fG.,  ;    is   taken,  anc   r's corded,   as    the 

abscissa  of   z-,  ^      Tae   uprjer  and   lox'er  bounds    of   s^   aro 

next   examined,    etCc       (Wg  h  s-o   assiMia   that 

R^z^    <  R^z^   <    '"    <  R    z      ^    <  R    z    ) -. 
Si—     e2~  —     e   n-"l  --     e  n 

Suppose  that  z-,  is  the  first  root  for  which  the 

upper,  Gt^,.  and  Icuer,  Q-r.    estir..ite3  of  the  f■b3f^issa  fail 

to  satisfy  the  inequality  <3  above.   The  progran  explained 

in  Section  2  then  finds  the  number  of  roots  to  the  left 

of  the  line 


z  =  •|(gK+Gj^)  +  i^j  -  00  < 


<   -T-   00  ; 


Gall   this   number    y)  ^      If  -^(gr^-'-G-^)    is    less    than 

G^  (v=K,  K-i-l,  .,. .  o,  ^  ),    it    replaces   G^,      •'■•^  2^'^K'^^K^    ^^    greater 
than  g    (v-\+l,\+2,  ■■  c  r.,n)    it    replaces   g^^  -      This   process 
is    repeated  \intil  g^,  and  G„  satisfy  Inequality  0  above c 
The   estimates    of   the    real  part   of   z^,-,    are    then   examined^ 
etc. 

In  this   way   the    abscissae   of   the   roots    z.,,,.  co.z 

J  1'         *    n 

are   determined    (usually  to  within  seven  sig.^ificant 
figures )  - 

The    ordinates    of   the   roots   are   now   found  by 
examining   the  polynomials    U    "(t)    and  V^    ' (t)    given  as 
the   real   and    imaginary  parts    of   f?   (R   z-^.+it)    - 
U^^^(t)    +   IV^^Mt)...      It    is    easy  to    see   that   u'-^^(t)    and 


V   ^(t)  have  Jlazy^   as  a  common  rsal  roo^.r   Unfortunately, 

U^  ^(t)  and  V^  ' (t)  will  have  as  common  roots  all  pairs 

of  roots  of  P  (z)  =  0  which  are  symmetrically  situated 

with  respect  to  the  line  z  =  B.   Zj^   ■•!-  it.   Therefore,  to 
^  e  K    • 

find  the  ordinate  of  the  root  z„,  the  greatest  common 
divisor  Q^^^  of  U^^^  and  V^^'  is  found.   The  real  roots 
of  Q   ^  are  foiind  by  determining  the  abscissae  of  its 


roots  and  recording  those  whose  imaginary  parts  vanish„ 

oef' 
(K) 


(K) 
(Q^  '  will  be  either  odd  or  even  if  ■■'j,(z)  has  real  coef- 


ficients and  it  is  possible  to  halve  the  degree  of  Q 

through  a  simple  change  of  variable,) 

Section  2-   Mathematical  Background 

In  order  to  locate  the  zeros  of  the  polynomial 

P  (z),  we  first  determine  a  large  circle  C  in  which  they 
^  max 


lie,   A  simple  bound  is  given  by  |z|  < 


1^1 

n 


il  +  1  +  10-^ 


where  ^-^(z)   =   a^z  +  "  •'   +   a,z  +  a^e   Application  is  then 
made  of  the  Cauchy  formula,  which  counts  the  number  of 
zeros  in  the  interior  of  a  simple  Jordan  curve,  in  a 
process  of  successive  subdivision  of  the  circumscribed 
square  around  C  ■,   In  this  way,  the  coordinates  of  the 
zeros  may  be  determined  theoretically  to  any  desired 
degree  of  accuracy. 

The  Cauchy  formula 


-b- 


(1)    N  =  -±-  jfdilog   ?„(z)] 


2711,'         ^ 


gives  the  number  of  roots  of  P  (z)  =  0  contained  within 
the  curve  B  (if  no  roots  are  on  B)r   It  would  be  possible 
in  principle  to  perform  a  nvmierical  integration  to  evalu- 
ate the  number  N,   But,  by  modifying  B  in  (1)  so  that  we 
obtain  a  formula  for  the  number  of  zeros  to  the  "left" 
of  any  straight  line,  it  becomes  possible  to  evaluate  the 
resulting  integral  by  simpler  algebraic  processes.   We 
now  give  the  details  of  the  above  program. 

Formula  (1)  counts  the  total  number  of  revolutions 
about  the  origin  made  by  the  radius  vector  drawn  from  the 
origin  to  the  point  P„(z)  as  z  varies  around  B,   Let  us 
assiome  B  consists  of  an  arc  of  a  circle  of  large  radius 
centered  at  the  origin,  together  with  its  chord.   Suppose 
that  the  chord  lies  on  a  fixed  straight  line  L  and  that 
the  circular  arc  lies  to  the  "left"  of  the  line  L.   If 
we  let  the  radius  of  the  circular  arc  increase,  the 
formula  (1)  will  yield  the  number  of  zeros,  q,  to  the  left 
of  the  line  L,   The  change  in  angle,  AJD,  of  P  (z)  =  Re"^ 
over  the  circular  arc  will  tend  to  nTi  as  the  radius  Increases 
since  for  z  =  re"''  and  r  large 

^n^^)  ^  V""  ■*•  ^n-l""""^  "■    '-'    ■•  ^1^  -^  ^0 
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n  "*      n 


+    ••  •   +  -) 


a  z 

n 


n 


^  a^z     =  a  r  e 


n  in© 


n 


n 


and    ^0  — >  Ti   as   r  — >  oo  • 


Therefore 


a   z 
n 


.n 


q  =  §  +  ~i-/^d   log    P„(z)    . 
27tl/ 


n 


The  principle   objective   of   the   present   section  is    to   show 
how   the  change    A®  in  the   angle    iD  as    z   varies   over  the 
line  L  may  be   expressed  algebraically.      To  this    end  we 
represent   the   points   on  the    line   L  by  z   =  y  +   t5,    where 
Y  and  5   are  fixed  complex  numbers   and   t   varies   from 
-  00  to  ■  +  00  .      Then 


Pj^(z)    =   Pj^{Y+t5)    =  U(t)    +   iV(t)    =  R(t)e 


ia)(t) 


(We  may  assume  that  the  polynomial  V(t)  is  of  degree 
>  the  degree  of  U(t),  otherwise  we  would  consider 
iP„{z)  which  has  the  same  zeros  as  ?^(^)»)      'I'he  number 
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of  revolutions  mads  by  !D(t)  can  be  found  by  appropriately 
counting  the  zeros  of  V(t),   Since  the  degree  of  V  is  not 
zero  and  is  >  the  degree  of  U,  we  know  that  !D(t)  as 
t  — >  +  00  tends  to  an  angle  which  is  not  0  or  7t,,   This 
means  that  if  we  count  the  number  of  ''positive"  crossings 
of  the  positive  and  negative  real  axis  made  by  the  radius 
vector  of  angle  !D(t)  and  subtract  the  niimber  of  "negative" 
crossings,  we  will  obtain  twice  the  number  of  revolutions 
made  by  the  radius  vector.   We  now  describe  how  the  count- 
ing may  be  performed  algebraically, 

A  "positive"  crossing  of  the  real  axis  in  the  (U, V) 
plane  is  one  for  which  V  increases  from  negative  to  positive 
values  if  U  is  positive;  while  if  U  is  negative  V  decreases 
from  positive  to  negative  values  as  t  increases.   The 
"negative"  crossings  are  those  in  which  V  has  the  opposite 
behavior. 

Now,  in  order  to  count  the  number  of  crossings; 
we  first  use  the  Euclidean  algorithm  for  finding  the 
greatest  common  divisor  of  V  and  U  in  the  form 

V  =  Uq-L  -  r^^ 

U  =  V-^q^   -    v^ 

« 
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where  degree  of  U  <  degree  of  V 
degree  of  r,  <  degree  of  U 
degree  of  r .  -,  <  degree  of  r. 

u  J 

If  we  now  consider  the  number  M(t)  of  changes  of 
sign  in  the  sequence  of  terms  (V, U,  r, ,  Tp, ,  ..r^,),  we  ob- 
serve that  M(t)  can  change  as  t  Increases  only  when  V 
changes  sign  since 

(1)  V  and  U  have  no  common  zero  (by  the  assumption 
that  the  line  L  was  zero  free) 

(2)  Whence  no  two  successive  terms  vanish  simultane- 
ously, 

(3)  Hence  the  nearest  neighbors  of  a  term  which 
vanishes  have  opposite  sign, 

(I4.)  Thence  the  only  way  for  M  to  change  is  by 

having  V  change  sign,  since  r„  cannot  vanish. 
Now  the  number  M(t)  decreases  for  a  positive  crossing  and 
Increases  for  a  negative  crossing. 
Hence 

M(-oo)  -  M(+oo  )  =  2  c  number  of  revolutions  of  (D  over  L 

whence 

_  n  ,  M(-oo  )  -  M(+co  ) 

_  -_  +  ^ 

2 

-10- 


^-2 


We  note   that   only  the   signs   and  the   degrees   of  the   lead- 
ing  coefficients   of    (V, U,  r-, ,  . ,  ,,r^)    are  needed  to  de- 
termine  the   numbers  M(-oo  )    and  M(+oo)^      (In   the   event 
that   P„(z)    has   real  coefficients,    and  the   line   L  is   a 
vertical   line,    one   of   the   functions   V, U  will  be   odd  and 
the   other  will  be   even;      The   sequence    (V,  U,  r, ,  , . .,  r~) 
will   therefore   alternately  be   odd  and  even  functions.      It 
is    therefore  possible   to   coinpute   M(-oo  )    -  M(+oo  )    from  a 
knowledge   only  of   the   signs   of   the    leading   coefficients 
of    (V, U, r,,  » , ,,  rr,) ,      For  example,    if   all   of   the   leading 
coefficients   have  positive   signs,    then  M(-oo  )    =  k  +  1 
while   M(+oo )   =   0,      In  other  words   if  the   leading   coef- 
ficients  of   two   neighboring  terras   have   the   same    sign,    they 
contribute   +  1   to  M(-oo )    and  0  to  M(+oo)    and  vice-versa 
if  the    adjacent   leading  coefficients   have   opposite   sign 
they  contribute   0  to  M(-oo  )    and  +   1   to  M(+co),      It   is 
therefore  possible  to   compute  M(-co)    -  M(+oo)   without 
storing   all   of  the   coefficients   of   all   of  the  terms 
(V,  U,  r, ,  , , ,,  r..)    in  the  process   of  performing   the  greatest 
common  divisor  algorithm.      That   is,    a  single    intermediate 
coiinter  can  be   used  to  store  K   (-oo  )    -  M   (+oo  ),    the   differ- 
ence   in  the  number  of   sign  changes   of   the  sequence 
(V, U,  r, ,  ,  .,,r    ),    while   the    term  r      ,    is   found  from  the 
knowledge   of   the    two  preceding   terms   r^_-]_   and  r^e 
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Aftc:'^  r.,,_.,    Ic.    foundj    r,.^-.    raaj  be   erased „      The   signs   of 
the    leading   coe.tTicient   of  ^^^...^    and   r     can  be    compared 
':-    find   the   numbei'  M.  ...  (-en)    »   M,^t(+0C')   vihich  reioresents 
the  difference   in  the   niimber  of  sign  changes    in  the 
sequence    (V.  U,.  r-, ,.  .-,r ->  r.  ,  r^. ,  -  )    etCo) 

Sect  ion  3--      Ma  chine   Aspects 

In  applying   the   method   discussed  in  the   preceeding 
sections    on  the   Univac    there   are    certain  difficulties, 
some    of  Khich   can   onl^r  partially  be   overcome o 

Consider  first   the   sequence    (V,  U,  r, ,  rp,  , ., ,  r^.) 
and  the   Euclidean  algorithm  which  leads    to   it 

^  "   '^2^-''l   -   ^2 
^1  "   "^3-2   "   ^3 


■K-2  "    ^iriv-1    "    "K 
K-l        %+1-K 


o 


The  degree  of  ^  is  knovn,  it  is  the  same  as  the  degree 
of  the  original  polynomial.   However,  in  order  to  success- 
fully carry  ov.t  the  .-^.Igoritlim^  one  must  be  able  to  determine 
correctly  the  degree  of  U  and  of  the  r'Sc   This  requires 
a  knowledge  of  what  terms  are  to  be  considered  as  being 
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zeroa      PurtheVj    it   ±z    ir-portaiit   to   know  \inat   is   to   be 
considered   as    zero   sc   as   to  Icnow  ;/hen  to   tsrminats   the 
algorith'T.,;      In  an  efforu   to    clrcvrnvent   this   difficulty 
an  artifice  was   devised-      Observe   first   that   each  member 
of   the    sequence    (V,  U,™  r, ;,  o -j  r;  r^.  ■    is    a  polynomial  which  con- 
tains    exclusively  odd  or   even  degree  terms.      We    then  assurre 
that    the  degree   of  U  is    one   less    than   the   degree   of  V,    If 
in  performing    the    necessary  division  this   assumption  leads 
to   an  exponent   overflow   or   a  zero  divisor  this    fact   is 
noted  on  the    supervisory  control   typevxriter    (SCT)    and   the 
division   is   performed  again  —   this   time   assuming  that    the 
degree  of  U   is    3    less    than   the   degree    of  V.,      This   process 
is   repeated  until  the   division  is    successfully  completedo 
An  entirely  similar  pi^ocedure    is    followed   throughout   the 
algorithm.      The   algoi'ithm  is    allowed  to    continue   until 
the    fin^l  terra  of  the    sequence  r„  is   a  constant    (degree 
zero  polynomial)  c      The    reason  for  this   is    clear,,      The    final 
term  of   the   sequence    should  be  the   greatesi:   comiaon  divisor 
(g,c.d    )    of  V  and  U,      This   goC«d.-,  xjIII  be   a   constant  pro- 
vided the   vertical   line   we   are   trying   does   not  pass    through 
a  root  exactly  and  prox^ided   there   are  no  roots    symmetrically 
placed   on  either  side    of  this    line,-      These   assumptions    are 
tacitly  made,, 
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3y- succQoi^ivo    applica'^icns    fs^^Qn.,.     (see   sscticn  1) 
W8    can  hope   to   find   the    o.bEoir.'sae   o:'    rho    ."?ocl.f%.,      'J?h3ro 
ara,    ho'^'eve??,    tv-o    things  x-:hich  m-'iy  prevent    tais.-.. 

(1/    VJe   may  try  a  vsi-tical    lino   and   got    the    rosuit 
that   there    arc;  'Ti  roct3    to   th?:;    Z.p-Ct   of    that   line;.      But 
from  previous    ocmputations   i-re  may  have   already   ''found'- 
more    than'']    roots   to    the    left   of   that   line        This    incon- 
sist-ncy   stops    the    Univac    and    the    type  .ut    '•  ^  N   error'- 
occurs,      Th^re   is    nothing  "lore   or-^    can  do   in   this    situation, 

(2)    The    exact   opposite   can   oco\.r,,      We  may  find   that 
there    are    ^   roots    to   the   left   of   a   vertical    line   when   in 
reality  there   are    loss   'i.han     '^    roots    there..      No   check  can 
be  made    for   this   error^    and,    in  general      if  it   occurs,    it 
will   lead   to  poor   accu'-.'acy  in  the    final   result:, 

Thete   two    errors    occur   only  in   the   neighborhood   ol 
a   root    and  particularly  in   the   neighborhood  of  a     multiple 
root.      In   an  effort   to   keep   the   oocii -.'rer.oe   of   th-^i;e   errors 
dc^n  to    a  mi.nimuj'tt  multiple   roots    a:r-e   oliiminated  at    the 
very  outset   of   the    computation  by  .■'.'actorlng  out   the   g,Cc,d,-. 
of   the    polynomial    and    its    derivative.       An  attempt   is    made 
to    compute   the   proper  multiplicity  of   each   root    (along  with 
an  error  estimate)    at    the   end   of    the   computation    (see   below). j 

Cnce   the    abscissas   hsve   been  found,    one    can  proceed 
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In  two   different,  wayt;    in  order   to   obtain   the   ordinatesc 
The   basic  me  shod    (Soctions    '.   and   2)    can  be   used  so  as    to 
yield  the   nu:;:ber  of  roots    aboj   or  below   any  given  hori- 
zontal  line  J      This,    however,    would   introduce   complox 
arithmetic   and   lead  to   extraneous    roots-      Therefore   a 
different  method  was    adopted-      Lot  1-    {?■)    be    the   polynomial 
that   wo    are    interested   in,    and   assvme   that;  we   already  know 
that  a   is   an  abscissa  and  that  we   vrish  to  find  a   correspond- 
ing  ordi'^ute   p^      I.eL   z   •=  a  ■»-   it     where   t    is    a   variable. 
Then  l^iz)   •■=  .".^^(c+it)      =  Uft)    -:-   iV't),        Then  clearly  we 
must  have   U{p)    -  V(6)    -   0.      Thj    oc-.r/crce   is    also    trv.3} 
I.e.,    if  U(p)    -  V(p)    =   0  for   real   (3   then  (3    is    an  ordinate 
corresponding   to   the    abscissa  a^      It    is    obvious    that    the 
real   roots    of   g(t),    the   goC^df    of   U(t}    and  V(t)j    are    the 
ordinates   that  vrc   desire.      We   notice   first   that   this    goC.d, 
g(t)    contains    only  odd   or   only  even  degree   terms ^      If   it 
contains    only  odd  degree   terms   then   a   factor   t   can  be   re- 
moved  and   thus   we   may  assume   that   g(t)    has    only  even  de- 
gree  terms  0      A   simple    change   of  variable  then  enables   us 
to  halve    the    dc^gree    of  the   g,_, Cod-      We  mus  :  then  find   the 
real   roots    of  this    '-reduced  g.c.dc"..    g-r,(t)»      If  the   degree 
o^  g^('t)    is    one,    the   results    can  be    rea":  right    off.      If 
the    degree   is    two,    a   simple    application  of   the    quadratic 
formula  gives    the   desired  result..      If  the    degree   is    greater 
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than  or  equal  to  three  the  typeout  "G„C.D,  DEG  6+"  occurs. 
Then  using  g  (t)  and  the  original  basic  method  (Sections 
1  and  2)  we  find  the  abscissas  of  g  (t).   An  e  -  test 
against  the  remainder  after  synthetic  division  is  then 
used  to  find  out  which  of  these  abscissas  are  real  roots. 
It  is  possible  (although  unlikely)  for  this  test  to  go 
astray.   If  this  happens  the  Univac  stops  and  the  type- 
out  "extra  roots"  occurs  at  a  later  time. 

The  main  difficulty  is  in  finding  the  g,c,d,  of 
V(t)  and  U(t),   The  Euclidean  algorithm  must  be  stopped 
at  precisely  the  right  step  —  not  too  soon  and  not  too 
late.   That  is  to  say  we  must  have  a  precise  way  of  telling 
when  a  division  yields  a  remainder  which  is  to  be  considered 
as  being  zero.   The  only  check  available  to  us  at  this  point 
is  that  we  know  the  g,c,d,  should  not  be  a  constant.   If  it 
is  calculated  to  be  a  constant  then  the  e  used  to  test  when 
a  division  yields  a  remainder  of  zero  is  multiplied  by  ten 
and  the  g,c,d,  is  re-calculated.   This  process  is  continued 
until  either  a  non-constant  g,c.d,  is  obtained  or  until  the 
original  e  has  been  multiplied  by  10  ,  at  which  time  the 
Univac  stops  and  the  typeout  "g,c,d,  const,"  occurs.  The 
g,c,d,  finally  obtained  still  may  not  be  the  correct  one. 
If  this  happens  either  the  Univac  will  stop  and  type 
out  "extra  roots"  at  a  later  tirn.e  or  the  error  will  manifest 
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itseir  ir.  bad  error  estimates  (see  below)  „ 

After  all  the  roots  of  the  reduced  (no  multiple 
roots)  polynomial  have  been  found,  we  try  to  obtain  an 
error  estimate  and  the  proper  multiplicity  for  each  root. 
We  proceed  as  follows:  Let  ■'?n(z)  be  the  original  un- 
reduced (with  multiple  roots  possible)  polynomial,  and 
let  a  +  ib  be  a  root.   Then  make  the  substitution 
z  =  a  +  it: 


P^(2)  =  ?^(a+it)  =  U(t)  +  iV(t). 

Divide  U(t)  synthetically  by  (t-b)  and  get  a  remainder 

r^*^^  and  a  quotient  q^   o   Divide  V(t)  synthetically  by 

(t-b)  and  get  a  remainder  rjL   and  a  quotient  q^  ', 

Then  ||  r,^  '  I  +  \r      '  u/a     where  a  is  the  leading  coefficient 

of  "^  (z)  is  the  error  estimate  used  for  the  root  a  +  ib, 
n 

Now  divide  both  q^°^  and  qj°^  synthetically  by  (t-b) 

and  get  remainders  r^  K    v^   ^    and  quotients  q^  ',    q^  ^ 

(11      (1) 
respectively.   If  both  r  '    and  V    '    are   less  than  some 
■^       "^  u        V 

computed  s,  then  a  +  ib  is  said  to  have  multiplicity 
twor   By  repealed  such  synthetic  divisions  we  may  hope 
to  compute  the  proper  multiplicity  of  each  rooto 

As  a  matter  of  actual  practice  things  do  not  always 
work  out  quite  so  nicely,  and  the  sum  total  of  the  multi- 
plicities may  not  add  up  to  the  degree  of  the  original 
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polynomial.   Particularly  in  the  case  of  fairly  high 
multiplicities  (say  >  3)  the  indicated  multiplicity  may 
not  be  correct  (generally  it  will  be  less  than  the  true 
one)c   A  multiplicity  of  zero  indicates  that  the  computed 
root  is  so  far  from  the  true  root  that  it  should  not  be 
considered  as  a  root  at  all. 

Many  of  the  machine  difficulties  mentioned  in  the 
preceeding  paragraphs  can  be  completely  or  partially 
avoided  by  scaling.   For  example  the  existing  routine 

Q   Q 

will  not  solve  properly  the  polynomial  P„(z)  =  10  z  + 

7  7      A  A 
10' z'  +  10  z  +••"+  lOz  +  1;  however,  no  trouble  is 

a  'y 

encountered  with  the  polynomial  ^^,(2)  =  z  +  z  +•••+  z  +  1, 
The  optimiom  condition  seems  to  be  one  In  which  all  the  coef- 
ficients are  of  the  same  order  of  magnitude,  and,  the  closer 
one  can  come  to  this  condition,  by  scaling  if  necessary,  the 
better  will  be  the  results  obtained^ 

The  time  taken  to  complete  a  problem  is  a  function 
of  the  degree  of  the  polynomial  and  of  the  geometry  of 
its  roots.   The  twentieth  roots  of  unity  take  slightly 


1,  Note:   If  the  polynomial  being  solved  is  given 
by  Pj^(z)  =  a^z^  +.»•+  a-j^z  +  aQ,  then  a  necessary  condition 
for  the  operation  of  the  routine  is  that  (max  |aj[|)^  di- 
vided by  (aj^)^~l  shall  not  cause  exponent  overflow  in 
floating  decimal. 
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over  t;i;o  hours;  the  tenth  roots  of  unity  takes  about  20 
minutes,^  a  tenth  degree  polynomial  with  all  real  roots 
takes  about  30  minutes,  etc. 
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Appendix  I  -  Operating  Instructions 

I,  Servos    and  Tapes. 

A,  Servo   5   -  Program  tape. 

B,  Servo  A  -  Blank  for   re-run  purposes, 

C.  Servo  6   -  Blank  for  edited  output    (optional) 

D.  Servo  7    -  Polynomial    (optional) 

II.  Breakpoints 

A.  Breakpoint  No.  Ic 

1)  Normal  Path  -  Polynomial  to  be  solved  read 
in  from  servo  7  and  servo  7  rewound.   The  tape 
must  be  prepared  in  the  following  manner:  The 
first  word  to  be  typed  must  be  OOOOOOOOOOnn 
where  n  is  the  degree  of  the  polynomial 
(01<nn<20),   The  next  word  must  be  the  constant 
term  of  the  polynomial,  the  next  word  the  de- 
gree one  term,  degree  two  term,  etc.,  up  to  the 
degree  n  term.   The  remainder  of  the  block  is 
immaterial.   All  coefficients  must  be  in  float- 
ing decimal  composed  form, 

2)  Force  Transfer  -  Computer  calls  for  the 
manual  type  in  of  the  polynomial  to  be  solved. 
In  response  to  the  type  out  "Ti  DEGREE"  type 
in  OOOOOOOOOOnn  where  n  is  the  degree  of  the 
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polynomial  (01<nn<20).   There  will  then  be  a 
series  of  type  outs  of  the  form  "TijjTERM" 
where  00<jj<nn,   In  response  to  each  of  these 
type  outSj  type  in  the  j-th  term  (floating 
decimal  composed)  of  the  polynomial, 

B,  Breakpoint  Noc  2o 

1)  Normal  Path  --  Final  edited  result  will  be 
written  out  at  20  pulses/inch  on  servo  6,  servo 
6  will  be  rewound,  and  the  number  of  blocks 
written  out  will  by  typed  out  on  SCT. 

2)  Force  Transfer  >-  Final  edited  result  will 
be  printed  on  the  SCT  (see  below  for  type- 
writer instructions) o 

C,  Breakpoint  NOo  9o 

Prom  time  to  time  servo  A  will  be  written  on 
and  then  rewoundc   To  re-run  from  that  point  (1) 
rewind  servo  5  ^^   necessary,  (2)  initial  read 
servo  5,  (3)  force  transfer  on  breakpoint  No.  9. 
III.  Typewriter  Instructions c 

A,  Always  keep  the  SCT  on  normal. 

B.  For  printing  final  edited  results  use  narrow 
paper  and  keep  both  left  and  right  hand  margins 
as  far  over  as  possible.   No  tabs  are  necessary. 
There  is  a  ^  at  the  beginning  of  the  editing  and 
3/' s  at  the  end. 

To  start  the  program  initial  read  servo  $., 

-21- 


Append5.x   II   •>  Sample    Problem 


Polynomial 


degre 
degre 
der"^e 
de^.e 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 
degre 


00 
01 
02 
03 

Ok 
05 
06 

07 
08 
09 
10 
11 
12 
13 
Ik 
15 
16 
17 
18 


term  -  +„100000000 
term  =  +^000000000 
term  =  +.000000000 
term  =  +o 000000000 
term  =  +,000000000 
term  =  +c 000000000 
term  •-  +,000000000 
term  =  +c 000000000 
term  =  +•»  000000000 
term  =  -o200000000 
term  =  +,00000000  0 
term  =  +„000000000 
term  -  +<, 000000000 
term  =  +.000000000 
term  =  +„ 000000000 
term  =  +^.000000000 
term  =  +.000000000 
term  =  +.000000000 
term  =  +.100000000 


(+01) 
(+00) 
(+00) 
(-^00) 
(-00) 
(+00) 

(+00! 

(+00) 
(-:00) 
(+01) 
(+00^ 
(+00) 
(  +  00) 
(+00) 
(+00) 
(  +  00) 
(  +  00) 

r+00) 

(+01) 


Roots 

real  part 

-c9396926l6  (+00) 

-o9396926l6  (+00) 

-,500000001+  (+00) 

-o5ooooooo[j.  (+00) 

+  ,17361+8179  (+00) 

+  .17361+8179  (+00) 

+  .7660I;14j.38  (+00) 

+  .7660i|ijii.38  (+00) 

+^999999996  (+00) 


imag.part 
+  ,3i+20201[j4 
-c3ii.20201[i.l+ 
+  .8660251^.18 
-c8660251|l8 
+  ,981+807732 
-.981+807732 
+.642787590 
-c  61^27  87  590 
+  ,0)0000000 


error 

mult 

(  +  C0) 

+  „8i;2  02  011+1+ 

(-07) 

02 

(^-00) 

+  ,81+20201^1+ 

(-07) 

02 

(  +  00) 

+  021+3683202 

(-06) 

02 

(  +  00) 

+  .21^3683202 

(-06) 

02 

(+00) 

+  ,-566855781 

("■07) 

02 

(  +  00) 

+  .■566  8557  81 

(-07) 

02 

(+00) 

+  -37953801^3 

(-06) 

02 

(  +  00) 

+  0  37  953  801+3 

(-06) 

02 

(+00) 

+  ,1+0000000  3 

(-07) 

02 

finished 
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